In-plane noncollinear exchange coupling mediated by helical edge states in Quantum 

Spin Hall system 
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We study the Ruderman-Kittel-Kasuya-Yoshida (RKKY) interaction mediated by helical edge 
states in quantum spin hall system. The helical edge states induce an in-plane noncollinear exchange 
coupling between two local spins, in contrast to the isotropic coupling induced in normal metal. The 
angle between the two local spins in the ground state depends on the Fermi level. This property 
may be used to control the angle of spins by tuning the electric gate. 



I. INTRODUCTION 

Recently the study of quantum spin hall (QSH) state 
has attracted much attcntionii^i^iii^ for its interesting 
topology and for its potential applications in the field 
of spintronics. 

The well-known examples of the topologically nontriv- 
ial states are the integer and fractional quantum hall 
states, where the quantization of the Hall conductance 
is protected by a topological invariant. The QSH insu- 
lator is a novel topologically insulating phase with time- 
reversal symmetry. A 2-dimensional QSH insulator has 
a charge excitation gap in the bulk and gapless helical 
edge states. Two states with opposite spin-polarization 
counter-propagate at a given edgeii^. Due to the time re- 
versal symmetry, these spin-filtered edge states are stable 
against weak interaction and disordes^. Hence, they can 
be viewed as effective spin and charge conducting chan- 
nels, and could be used to construct promising spintronic 
devices with low power consumption. 

RKKY interaction is an effective interaction between 
two local spins mediated by conduction electrons^. It 
plays an important role in many fields of solid state 
physics, e.g. giant magnetoresistance (GMR)^, dilute 
magnetic semiconductor (DMS)^. More recently, peo- 
ple propose that the controllable RKKY interaction can 
be used to manipulate the quantum states of the lo- 
cal spins, a crucial point for the spintronics and quan- 
tum computingi^iii. Thereafter, the RKKY interac- 
tions in different spintronic materials, e.g. spin-orbital 
systemi^ii^iii^ and graphen o^^d^d'^d^ , have been investi- 
gated carefully in order to facilitate the further develop- 
ment of the spintronic devices. Compared with the nor- 
mal cases, the effective interactions between local spins 
in these systems exhibit rather different properties. The 
RKKY interaction in the spin-orbital system becomes a 
twisted exchange coupling since it is a spin-dependent 
systemi^ii^iii^. The spin polarizations of the two local 
spins are no longer coUinear in this case. As for the 
graphene, due to its special electronic band dispersion, 
it is found that the RKKY interaction is ferromagnetic 



for local spins within equivalent sublattices but antifer- 
romagnetic for opposite sublattices when the Fermi level 
is near the Dirac poin t^^d^ . 

Being an novel spintronic material, it is an quite in- 
triguing and practical problem that what the RKKY in- 
teraction in the QSH insulator is. Since the helical edge 
states are the only conducting channels in this system, 
the problem becomes that what the exchange interaction 
mediated by the helical edge states is. 

In this paper, we investigate the RKKY interaction 
mediated by the helical edge states in the QSH insula- 
tor. Our theoretical analyzing is mainly based on the 
simplified model of the helical edge state^i^i^ , in which 
the spin of the carriers is assumed to be parallel or an- 
tiparallel along the z axis as shown in Fig. 1. Actually, 
it is proposed that the QSH insulator could be realized 
in various kinds of systems, e.g. graphene^, quantum 
well^i^, or semiconductor materials with special strain 
gradient^i^. The descriptions of the helical edge state in 
different systems may be different. However, the model 
we used is the simplest and most basic one. It grasps the 
primary characteristic of the helical states, i.e. the corre- 
lation between the spin polarization and propagation of 
the carrier, and also is the exact expression of the heli- 
cal edge states in the semiconductor system with special 
strain gradient. Hence, we believe that this model is a 
good starting point for analyzing of the RKKY interac- 
tion mediated by the helical edge state. 

We find that the helical property and the linear dis- 
persion of the edge conducting electrons will lead to an 
in-plane and noncollinear exchange coupling between two 
local spins along the edge. In the asymptotic limit, this 
interaction has a simple expression. We can see that the 
angle between the two spins can be controlled by adjust- 
ing the the Fermi energy of the system. When the Fermi 
surface is near the Dirac point, i.e. the Fermi energy is 
near zero, the effective interaction becomes a constant 
antiferromagnetic exchange coupling. Actually, due to 
the helicity, pure spin current can be achieved in these 
edge states. Hence, this effective coupling can also be 
viewed as RKKY interaction mediated by pure spin cur- 
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FIG. 1: Schematic of the RKKY interaction between two local 
spins mediated by the helical edge states in QSH insulator. 
The QSH insulator is assumed to be a finite strip geometry 
which is infinite along the x direction. Two local spins Si and 
S2 are located along one edge of the QSH insulator. The solid 
(dashed) line represent the up spin right movers (down spin 
left movers) of the helical edge states on this edge. 



rent. 

This paper is organized as follows. In Sec. II, the 
RKKY interaction mediated by the helical edge states 
of QSH insulator is derived. In Sec. Ill, we discuss the 
special properties of this exchange interaction. Finally, a 
brief summary will be given in Sec. IV. 



II. MODEL AND FORMALISM 

We consider the simplest model of helical edge states 
which has been successfully used in the study of the tun- 
neling properties of the helical edge states in the QSH 
insulator"^ '^^i^° . The schematic is shown in Fig. 1. The 
helicity correlates the spin polarization with the propa- 
gation. Here, we assume that the right (left) movers ipni 
(ipLi) carry spin up(down). In the noninteracting case, 
the linearized Hamiltonian is 



VF / dx(-0j^i9a:?AflT " '^tl^'^^^Ll) (1) 



Ho 



where is the Fermi velocity. 

The localized spins are magnetic impurities and are 
denoted by Si and 82- Normally the Kondo coupling 
between the local moments and conducting electrons are 

Hi^-'^ ^ (Jap-Si I dx^p:^{x)S{x - Xi)iJi3{x) (2) 

^ »=1.2 

where J is the coupling constant, (Tq^ is the spin oper- 
ator of the conducting electrons, a and /3 are the spin 
indices, i = 1, 2 is the index of the local spins since in or- 
der to study the RKKY interaction, we need to consider 
two local spins . For the helical edge states, the only 
difference is the spin polarization is correlated with the 
propagation. Hence in the Kondo coupling expression 
of the hehcal edge states, V't(^) ii'li^)) means ^Jr^{x) 



If the coupling J is small. Hi can be treated as a per- 
turbation on Hq. The RKKY interaction between two 
local spins Si and S2 can be calculated from the second 
order perturbation thcoryi^iii^i^^ 



E 



RK 



Im / dioTr[{Si ■ a) 



^ J-00 ^ (3) 
xG^Ri2;uj){S2-(T)G'-{~Ri2-,oj)] 

where Ep is the Fermi energy, R12 = xi — X2 is the dis- 
tance between the two local spins and Tr means the trace 
over the spin degree of freedom of conduction electrons. 
Therefore, the study of RKKY interaction has been re- 
duced to the calculation of the retarded Green's function 
of the helical edge states. 

The definition of retarded Green's function is 

Glp{xt,x't') = ~i6{t-t'){{i^^{xt),^+{x't')}) (4) 

where a,l3 — {R |,L |} are the helical spin indices. 
The calculation of Green's function is quite straightfor- 
ward. Because that there is no coupling between different 
spins in the noninteracting Hamiltonian ([1]) , the nonzero 



Green's functions are G 



{x,x') and G^'^^^^ix^x') 



Take G^j^^j^ (a:, a;') for example. 



G^j,^j,^{x,x'-Lo)^ / die-*GJ^T^T(xi,x'i') 



= ^ I dk ■ 

27r ./ uu — vpk + IT] 



(5) 



where A is the cutoff of the momentum and CA(|fc|) is the 
cutting off function. Since our model is a low energy ap- 
proximation, a cutoff of the momentum is necessary and 
according to the discussion in the study of Graphene^^, 
a sharp cutoff is not suitable here. Therefore, we use a 
smooth cutting function 



CA(|fc|) = e-V. 
With this cutting off function, we get 



(6) 



G 



mm{x, X -jUj) — ■ e 

Vp 



■ e "-f 



9{x-x') (7) 



where 9{x — x') is a step function. The retarded Green's 
function of the left movers is 



G'LiLi{x,x']uj) = e" 

Up 



.e'i^("'-")0(x'-s) (8) 



Here, e is the decay factor of the retarded Green's 
function. We have to emphasize that the Green's func- 
tions with different smooth cutting functions will have 
similar form. The only difference is the decay factor. 

Then substituting the Green's functions into Eq. ([3]), 
we will get the final expression of the exchange interac- 
tion between two local spins. Without loss of generality. 
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we set i?i2 = xi ~ X2 > and then 



Erk = Im / ^^[(5"! • cr)iTG^T-RT(-^i2) 

^ J —OG 

= Ini / dLj{SlxS2x + SlyS2y + iSlyS2x 

71" J-oo 

- «S'la;<S'2j/)G5j|j^|(i?i2)G£|^|(-i?12) 

(9) 

It is clear that the range function for terms of SixS2x 
and SiyS2y is 

Fi(i?i2) = -^Im da;G5s^^^(i?i2)GL^i^(-i?i2) 



(10) 



The range function for term of SiyS2x is 
t2 r^F 



t2 i-^f 

F2(i?i2) = Im[* / <kjG^j,^j,^{R,2)Gl^Li{~Ri2)]- 

(11) 

And that for term of SixS2y is 



^3(-Rl2) — —F2{Rl2) 



(12) 



Hence, we only need to consider the range functions 
Fi[Ri2) andF2(i?i2). 

We see that all the range functions are related to the 
kernel function 



K{R)^ dwG'-R^R^{R)Gl^L^{~R). (13) 

J —OQ 

It depends on the Fermi energy 



K{R) = 



2VF(\+iR) 



2VF L . 



Finally, for cases E'f < 0, we get the range functions 



(14) 



F,{R) ^ - 



sin(^^^) 2EpR 



F2iR) 



2nvFij2+R^) A 



vf 



cos{ 



2EfR\ 



27rvj.(i + i?2)^ A 



(15) 

Rsin{ — — — )] 



Vf 



And for cases Ef > 0, 



2(-Ef) 



Fi{R) = ^ 



\2t-EF)R] 



-{- 



2nvfi^+R^)' A 
2{-Ef)R, 



(16) 



(17) 



Rcos 



Vf 



-]} 



F2{R) 



2t:vf{ ^ 
cos( 



A2 

2EfR 



A 



2 _2«F 
Rsm( )J} 



(18) 



III. DISCUSSION 

In the continuum limit A|a; — x'| = oo, we will get a 
rather simple expression 

27rui.|i?|^ ^ ^ (19) 

-sina ■ {SiyS2x - SixS2y)] 

where R = x — x' is the distance between the two local 
spms and a = -'-^ — . 

If we consider the two local spins and S2 as classical 
spins, Eq. p9|) can be transformed into 



Erk — 



2TrvF\R\ 



sin9isin92 ■ cos[{(j)i — 4>2) + ct\ (20) 



Here, (M, 6^1,01) and (M, 02,02) are the spherical coor- 
dinates of the spin vector Si and 5*2. We can see that 
only if 6*1 = 62 = it/ 2 and cos[{(f>i — (/)2) + a] — — 1, 
the system will have its lowest energy. In this case, 
61 = 02 = tt/2 means that Si and S2 are in-plane. And 
cos [(01 — ^2) + ct] = ^1 shows that the exchange interac- 
tion is noncoUinear. It means that the effective exchange 
interaction mediated by the helical edge states is an in- 
plane and noncoUinear coupling. Actually, it is easy to 
see from the expressions that in-plane and noncoUinear 
are the general characteristics of the helical edge states 
mediated exchange interaction which does not depend on 
the continuum limit. 

In the continuum limit, the angle between the local 
spins is determined by a = 2\r\Ef ^ ^ ^j^^ Fermi en- 
ergy Ef and the distance R. Especially when the Fermi 
energy Ef = 0, i.e. the Fermi level is around the Dirac 
point, a = and 0i — 02 = tt. It means that the exchange 
coupling becomes a constant and is always antiferromag- 
netic. Here, the Dirac point is the crossing of the bands 
of right and left movers. However, in general cases, we 
do not have a simple formula about the angle between 
the local spins. It should be determined though concrete 
numerical calculation. 

The special characteristics of the helical edge states 
mediated exchange interaction result from the interplay 
between the helicity and the linear band dispersion. As 
shown in former studies, without helicity, if the system 
is spin- independent, the coupling will have similar form 
F{R)Si ■S2 and the only difference is just the range func- 
tion F{R). In our case, the helicity makes the system 
spin-dependent: though the matrix of Green's function 
is stiU diagonal but G5j|^|(a;, a;'; w) 7^ G2|^|(a;, a;'; w). It 
is the main reason of the noncoUinear behavior of the 
exchange coupling. In addition to the helicity, the linear 
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band dispersion induces opposite step functions into the 
Green's functions as shown in Eq. ([7]) and ([5]). Actu- 
ally, the in-plane characteristic of the exchange interac- 
tion mainly results from these step functions. 

IV. CONCLUSION 

In summary, based on the simplified model of the he- 
lical edge states, we have investigated the helical edge 
states that mediate RKKY interaction between local 
spins in the QSH system. Since the conducting elec- 
trons in the helical edge states are the only carriers in 
the QSH insulator, this exchange interaction is probably 
the only possible mechanism of the exchange coupling 
between local spins in such a promising spintronic sys- 
tem. Furthermore, due to the helicity, i.e. the correlation 
between the spin polarization and propagation, it is be- 
lieved that pure spin current can be realized in this edge 
state. It means that this exchange interaction is actually 
an exchange interaction mediated by pure spin current. 
Hence, this problem is not only of fundamental interest 
but also useful for the future development of spin-based 
devices in such systems. 

We analyze the simplest theoretical model of the helical 



edge states and concentrate on the effects of the helicity 
and its linear band dispersion. We find that the RKKY 
interaction mediated by this helical edge state in such a 
system is in-plane and noncoUinear, which is extremely 
different from the exchange interaction in other systems. 
The angle between the local spins depends on the Fermi 
energy of the system. Therefor, this effective interaction 
offers a possible way to control such angle through ad- 
justing the fermi level by a gate. In the continuum limit, 
a simple expression of this exchange interaction can be 
achieved. Especially, when the Fermi level is around the 
Dirac point, the exchange coupling becomes a constant 
antiferromagnetic one. We also point out that these pe- 
culiar properties result from the interplay between the 
helicity and the linear band dispersion. 

However, our study is based on the noninteracting low 
energy approximation model of the helical edge states. 
For a concrete material or experimental setup, many 
more practical factors may need to be included. We be- 
lieve that our analysis is a good starting point for the 
further investigation. 
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